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We report the design and measurement of Superconducting/normal/superconducting (SNS) prox- 
imity DC squids in the long junction limit, i.e. superconducting loops interrupted by two normal 
metal wires roughly a micrometer long. Thanks to the clean interface between the metals, at low 
temperature a large supercurrent flows through the device. The dc squid-like geometry leads to 
an almost complete periodic modulation of the critical current through the device by a magnetic 
flux, with a flux periodicity of a flux quantum h/2e through the SNS loop. In addition, we examine 
the entire field dependence, notably the low and high field dependence of the maximum switching 
current. In contrast with the well-known Fraunhoffer-type oscillations typical of short wide junc- 
tions, we find a monotonous gaussian extinction of the critical current at high field. As shown in 
[l5|]. this monotonous dependence is typical of long and narrow diffusive junctions. We also find 
in some cases a puzzling reentrance at low field. In contrast, the temperature dependence of the 
critical current is well described by the proximity effect theory, as found by Dubos et al. [ljf| on SNS 
wires in the long junction limit. The switching current distributions and hysteretic IV curves also 
suggest interesting dynamics of long SNS junctions with an important role played by the diffusion 
time across the junction. 
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Introduction 



The proximity effect, or the penetration of super- 
conducting correlations in a normal (non superconduct- 
ing) conductor due to the proximity to a superconduct- 
ing one, has been shown to give the normal conductor 
superconducting-like properties over a length that can 
be considerable at low temperatures. This is because at 
low temperature both the phase coherence length and 
the thermal length \JTiD /ksT can extend beyond the 
length of the normal metal. This quantum coherence 
of a mesoscopic conductor then ensures that Andreev 
pairs, electron-hole time-reversed pairs formed at the 
SN interface, can propagate over distances of the or- 
der of the phase coherence length in the normal metal. 
Such a length is notably greater than the superconduct- 
ing coherence length y/HDjA. The normal conductors 
have ranged from simple noble metals to semi-conducting 
planes or wires, to long or short molecules (carbon nan- 
otubes , endohedral fullerenes, DNA ) down to individual 
atoms in break junction geometries [H, 0, 

One of the most remarkable consequence of the prox- 
imity effect is the propensity of a normal conductor well 
connected to two superconductors with different phases 
to carry a supercurrent, thereby demonstrating that An- 
dreev pairs carry superconducting correlations along with 
information about the macroscopic phase of the super- 
conductor at the boundary. Depending on the length of 
the normal part, the maximal supercurrent (called crit- 
ical current I c ) that can be carried is determined cither 
by the superconducting gap of the superconductor, or 
by the characteristic energy of the normal metal, the 
so called Thouless Energy, which in a diffusive conduc- 
tor of diffusion constant D is given by Eth = hD/L 2 . 
Indeed the critical current has been shown fijl to be 



given by ei?Ar/ c (T = 0) = 2.07A in the case of short 
junctions [L < £5, or equivalently E^h > A), and 
ei?Ar/ c (T = 0) = 10.82ETh in the opposite limit of long 
junctions. More generally, the remarkable point is that 
in the short junction limit the superconductor determine 
the properties of the proximity system, and in the long 
junction limit it is the normal metal (via its length, diffu- 
sion time, coherence length) which determines the prox- 
imity effect. 

The interplay of the mesoscopic phase coherence in the 
normal metal and the superconducting coherence due to 
the nature of the superconducting state, have led to the 
fabrication of many interferometer devices. The experi- 
ments that launched the development of the mesoscopic 
proximity effect mainly consisted of a normal metal loop, 
with nearby superconducting islands ([!, 0]). The resis- 
tance of these devices was modulated by the magnetic 
field one or two orders of magnitude more than the uni- 
versal conductance fluctuations in normal rings with no 
superconductor. But no supercurrent was measured be- 
cause normal leads connected the mesoscopic devices. In- 
terferometers in which a supercurrent was obtained and 
modulated are much more recent, and in most cases do 
not use a noble metal as the normal part: the non su- 
perconducting part actually is made of a long molecule 
such as carbon nanotubes Q , or semiconducting wires 
In some cases tt junctions can be produced, either by the 
use of a ferromagnet between superconductors or by the 
action of a gate electrode which modulates the popula- 
tion of the normal part. 

Interestingly, although many mesoscopic interferome- 
ters have been devised with a normal metal between sev- 
eral superconductors ( 0, 0] ) , no equivalent to a dc-squid 
(a superconducting loop with two insulating barriers or 
weak links) has yet been fabricated in the mesoscopic 



regime [8]. In this article, we describe several proximity 
dc-squids, i.e. superconducting loops with two normal 
metal bridges, ranging from 0.75 to 1.9 microns long. 
We find that the IV curves are hysteretic, a phenomenon 
that we attempt to explain by the intrinsic dynamics of 
the diffusive normal metal. We find that the temperature 
dependence of the critical current of each long junction 
is in accordance to the theory of the proximity effect, 
much like in the work of Dubos et al. [16}. The field 
dependence is more original: In the proximity dcsquid 
configuration the critical current is highly modulated by 
the magnetic field. In both the single SNS junction and 
the dcsquid configuration, at large perpendicular mag- 
netic field, the extinction of critical current with field 
is monotonous and does not show the interference pat- 
tern previously observed in wider junctions @, 0, Il0| . 
We present two theoretical models which attempt to ex- 
plain this behaviour by including the diffusive nature of 
the normal metal. A surprising reentrance of the critical 
current at low field is also observed in some samples, and 
remains unexplained. 



B. Sample fabrication and measurement apparatus 

As a normal metal we used 99.9999% pure gold, with 
a content in magnetic impurities (Fe) smaller than 0.1 
part per million. Weak localization measurements in a 
simple wire made of this same material (at a different 
time) determined a phase coherence length of the order 
of 10 //m below 50 mK [20j. Two different superconduc- 
tors were used, either aluminum or niobium (see fig. [1] 
and [2}. The proximity wires and dc squids made from Al 
and Au were obtained by conventional double angle ther- 
mal evaporation through a suspended mask, in a vacuum 
of 10~ 6 mbar. The process with niobium as a supercon- 
ductor is more involved. We start with a bilayer of 50 
nm-thick thermally evaporated gold, covered by a 200 
nm-thick layer of sputtered niobium, and cover the bi- 
layer with a 70 nm-thick aluminum layer. With e-beam 
lithography we define and wet-etch open a window in the 
Al mask through which we locally etch away the Nb us- 
ing an SF6 plasma. In a realignment step we create a 
ring or wire shaped Al mask around which the bilayer is 
removed by an SF6 reactive ion etch followed by an ar- 
gon ion beam etch. Finally the Al mask is removed and 
the samples are tested before cool-down. The sample pa- 
rameters, obtained via SEM visualization or transport 
measurement, are gathered in table [U 

The length of the normal metal, which varies between 
0.75 and 1.9 fim, is much greater than the superconduct- 
ing coherence length £g = hD/A, where A is the BCS 
superconducting energy gap (see ratio Exh/ A in the ta- 
ble). Thus the devices are in the long junction regime, 
in which the properties of the normal metal are known 
to determine the proximity effect, rather than the prop- 
erties of the superconductor f [i~2|). 

The two superconductors used, Nb and Al, differ by 
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FIG. 1: SEM pictures of the two Nb/Au/Nb SNSdcsquids. 



FIG. 2: SEM pictures of Al/Au/Al SNS wires and SNS dc 
squids: from top to bottom, and left to right: wires Alwl250 
and Alwl300, SNS dc squids Alsql500, Alsq900 and Alsql900. 
Scale bar is 2 p,m. 

their transition temperature. The fabrication procedure 
also produces different S/N interface conditions: In the 
Nb/Au samples, the superconductor is a homogeneous 
bilayer of Nb/Au and the normal metal is a bare Au 
region. In contrast, in the angle-evaporated samples of 
Al/Au, the bare Au wires contact a region where Au 
and Al overlap (see SEM image). In addition, the Al 
loop is at some points, (depending on the loop shape, see 
picture) covered by Au. Therefore the superconductor in 
this configuration is less homogeneous. 

The geometry of the SNS long: junction is similar to 
previously measured ones 

(EE Hi)- The dcSNSquid ge- 
ometry has to our knowledge not yet been reported in 
mesoscopic samples. A dcSNSquid with sub millimeter- 
wide junctions was however produced by Clarke and Pa- 
terson in 1971 Q, in the aim of realizing an amplifier via 
a tuning of the dc squid asymetry. This will be detailed 
further on. The SNS dcsquid consists of a superconduct- 
ing loop interrupted by two normal metal wires. In our 
setup, superconducting leads in a four wire configuration 
connect the loop. In comparison, the usual dc squid is 
formed by a superconducting loop interrupted by insu- 
lating tunnel barriers. A dc supercurrent flows through 
such a dcsquid, and its maximum value, the switching 
current, is known to be modulated periodically by a mag- 
netic field, with a period of one flux quantum $o = ^/2e 
through the loop area. 

The measurements were conducted via room temper- 
ature 7T filters, home made cryogenic lossy coaxial cables 
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Nb750 


Nbl200 


Alw900 


Alwl300 


Alwl250 


Alsq900 


Alsql500 


Alsql900 


geometry 


SNS dc squid 


SNS dc squid 


SNS wire 


SNS wire 


SNS wire 


SNS dc squid 


SNS dc squid 


SNS dc squid 


Ljv(/ura) 


0.75 


1.2 


0.9 


1.3 


1.25 


0.9 


1.5 


1.9 


W(pm) 


0.4 


0.4 


0.125 


0.125 


0.125 


0.13 


0.15 


0.2 


A 


16 K 


16 K 


2.8 K 


2.8 K 


2.8 K 


2.8 K 


2.8 K 


2.8 K 


E Th from high T fit 


140 mK 


50 mK 


33 mK 


36 mK 


44 mK 


47 mK 


30 mK 


26 mK 


A/E Th 


114 


320 


85 


78 


64 


60 


93 


108 


£s = L N ^/E T h/A(nm) 


70 


67 


100 


150 


156 


116 


155 


180 


R N from dV/dl (fi) 


0.4 


0.7 


5 


5 


6 


2.6 K 


3.3 


3 


/r aa; (M) 


330 


68 


13 


3.7 


4.3 


17 


7.6 


4.6 


/r n (M) 


170 


26 


- 


- 


- 


9.2 


1.8 


0.6 




75 


16 


1 


0.6 


0.6 


1.3 


1 


0.45 


b 


11 


11 


22 


6 


7 


11 


10 


6 




4.4 


4.25 


13 


6.2 


7.2 


13 


7.6 


10.2 


\/2b 


4.7 


4.7 


6.6 


3.5 


3.8 


4.7 


4.5 


3.5 


T* (mK) 


100 


40 


35 


20 


25 


35 


21 


14 



TABLE I: Characteristics of the different samples. Ljv and W are respectively the length and width of the normal metal, which 
was 50 nm-thick high purity (99.9999%) gold for all samples. Eth = TiD/L 2 N is the Thouless energy, with D the diffusion 
constant. Note that £s = \JfiD j A is the superconducting coherence length determined from the gap of the superconductor 
but using the diffusion constant in the normal wire, so that it is evaluated via £s = L N ^E Th /A. The Al SNS wires (labeled 
Alwx, with x the length of the normal wire in nm) were produced on the same chip, and therefore should have the same 
material constants. This is also true for the Al SNS dc squids (labeled Alsqx), and for the Nb SNS dc squids (labeled Nbx). 
The superconducting gap A is deduced from the transition temperature of the superconducting contacts, using A = 1.76T C 



leading to the dilution stage, and 150 pF parallel capac- 
itors on the sample holder. The samples were current- 
biased with a dc current up to 400/xA and a small ac 
current (a few hundred nA) at a few hundred Hertz for 
differential resistance measurements. 



C. Transition temperatures 

The resistance of two Au/Nb dc squids is plotted in 
Fig. [3] as they are cooled to low temperature. The first 
resistance drop corresponds to the transition to a super- 
conducting state of the Nb/Au bilayer, around 9.5 K. Be- 
low the transition temperature, the resistance varies by 
less than 10% over a large temperature range. The final 
resistance drop to a zero resistance state, the sign that 
the proximity effect develops through the entire normal 
metal so that a supercurrent flows through the normal 
metal, occurs at a lower temperature, roughly between 
20 and 40 Etk depending on the samples. This tem- 
perature corresponds to the temperature at which the 
normal wire becomes phase coherent, and the proximity 
effect can set in. It is thus expected that it should cor- 
respond to min(Lr,L0) = L/2, i.e. ksT = AExh m the 
limit where Lt is smaller than L^. The fact that the 
transition temperature is 5 to 10 times larger than this 
limit recalls the factors 3 and 10 in front of the Thouless 
energy in the expressions of the induced minigap or the 
critical current decay energy (see related paragraphs). 
The exact 'proximity effect transition temperature' has 
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FIG. 3: Resistance of the two Nb/Au proximity dc squids 
versus temperature during cool-down. The configuration is 
a 4- wire measurement. The first transition around 9 K is 
the superconducting transition of Niobium. The transition at 
lower temperature (around 5 K for the shortest wire and 2 
K for the longest one) is the transition of the gold wire from 
normal to a proximity-induced superconducting state. 



to our knowledge not been calculated theoretically. 

The behavior of the three Al/Au SNS wires, plotted in 
Fig. Q]is similar. 
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FIG. 4: Resistance of three Al/Au/Al wires as a function of 
temperature during cool-down. Sample Alwl250 is measured 
in a four wire configuration, Alwl300 and Alw900 were mea- 
sured in a three wire configuration and a constant resistance 
of respectively 79 and 72 Q, corresponding to the resistance 
of one lead, has been subtracted. The double transition at 1.6 
K and 1 K for samples Alw900 and Alwl300 is due to their 
geometry, in which part of the superconducting wire is made 
of pure Al with a critical temperature of 1.6 K, and part is 
made of an Al/Au bilayer, with a critical temperature of 1 K. 



FIG. 6: Differential resistance of the shortest Au/Nb dcsquid 
(Nbsq750), at maximum and minimum switching current. 



100 
> 



50 





1.2 k y 


> 




> 






/ KmA) 



0.06 K 



50 




I(jiA) 



400 



FIG. 5: Typical V(I) curves of SNS dc squids at two magnetic 
fields, one in which the interference between both branches 
is constructive (3>/$o = 2n7r) and the other destructive 
(<&/<&o = (2n + 1)tt), showing a maximum and minimum 
switching current. The hysteresis is clearly visible and charac- 
terized by a large ratio (roughly 3) between switching current 
Is and retrapping current I r . The sample measured here is 
the shortest Au/Nb dc squid Nbsq750. The V(I) curves for 
SNS wires are similar. Inset: IV curve of the same sample 
at 1.2 K,when there is no hysteresis. The shape of the curve 
corresponds to the square root dependence expected in the 
overdamped regime, Q < 1 (continuous line). 



ence between both branches is constructive ($/$o = or 
multiple of 2n) and the other such that the interference is 
destructive ($/$o = (2n + l)7r). Single SNS junctions in 
the long limit have similar IV curves (not shown) . As the 
current is increased, a supercurrent flows (zero voltage 
drop), up to a value called the switching current (noted 
I s ), above which a voltage drop appears, and the junc- 
tion is said to have switched to the dissipative branch. 
As the current is decreased, the switch from the dissipa- 
tive branch to the superconducting one occurs at a cur- 
rent, called the retrapping current (noted I r ), which is 
much smaller than the switching current: the I V curves 
are strongly hysteretic, as found in previous experiments 
([l|| HH, [2y|). The value of the retrapping current does 
not depend on how far we sweep the DC current, nor does 
it depend on the temperature, in stark contrast with the 
switching current, as described in the next section. 

Equivalently, the differential resistance as a function 
of DC current is shown in fig. [6] for an SNS DC squid. 
The resistance in the normal state varies with current by 
a few percent in the Nb samples, and much more in the 
Al samples, probably because of inhomogeneities in the 
Al and Al/Au superconducting leads, which can generate 
phase slips. Some features on the down curve before re- 
trapping (strongly fluctuating differential resistance just 
before retrapping, and resistance bump at 140 /i^4)are 
still not fully understood, but should be explained by the 
dynamics of an SNS junction out of equilibrium [l6l. |28| . 



D. IV curves and differential resistance 

1. Experimental curves at low temperature 

Typical V(I) curves are shown in fig. [5] for a SNS dc 
squid at two magnetic fields, one in which the interfer- 



In the following paragraphs, we address the questions 
of the value of the switching current, which is 10 to 20% 
smaller than the critical current of the junction, the de- 
gree of hysteresis in the IV curves, and the fluctuation 
in the switching current for a given junction, field and 
temperature. 
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2. Dynamics of the SNS junction: the RSJ model adapted 
to SNS junctions 

The IV curves are well understood in SIS tunnel junc- 
tions, which are modeled by a Josephson element in par- 
allel with a capacitance and a resistive element (the so- 
called resistively shunted junction model, see for instance 
ref. [HI). The elements of the circuit lead to a dynamical 
equation for the superconducting phase of the junction (p 
given by: 

Cp+ — (p + u 2 p {sin{ip) - — ) = (1) 



Here io v = l/yf LrC where Lk is the kinetic induc- 
tance given by Lk = 4>q/2ttI c , and tr = RC controls 
the charge relaxation through the junction. The state 
of the junction is represented by a particle whose po- 
sition is the phase evolving in a tilted oscillating po- 
tential U(ip) = —Ej(cos(p + I/I c ip). The slope is the 
biasing current and the height of the barrier is related 
to the Josephson energy Ej — (h/e)I c . The time scale 
tr = RC represents friction for the particle. In the non 
dissipative state the particle is trapped in a local mini- 
mum, and oscillates at the plasma frequency uj p /2tt. In- 
creasing the current increases the slope of the potential 
and decreases the barrier height. Eventually the par- 
ticle escapes at I s < I c . The escape being a stochas- 
tic process, I s is characterized by a certain dispersion 
which is related to the thermal activation or the quan- 
tum tunneling of the particle through the barrier. The 
hysteresis is also understood within this model. On the 
dissipative branch, when decreasing the current below 
I s , the current consists in the externally imposed dc cur- 
rent plus an ac contribution due to the Josephson rela- 
tion ip = eV/h. At high dc voltages and currents, the ac 
part of the current, oscillating at the Josephson frequency 
u>j = h/2eV, flows entirely through the capacitance. The 
dc part of the current is dissipated in the resistance, and 
the IV curve is linear. At lower voltage, corresponding to 
Josephson frequencies approaching the circuit's resonant 
frequency u) p , the ac current starts to dissipate in the 
resistance, giving an additional contribution to the IV 
curve. The condition uj p — uij thus corresponds to the 
onset of non linearities in the IV curve and can be chosen 
as a criterion for the retrapping condition [27| , leading to 
I s /I r ss RCuip = Q. This explains that the retrapping 
current can be much smaller than the switching current 
if tr > l/ujp or Q > 1. 

a. Damping and hysteretic IV curves However this 
model is not directly applicable to the experiments on 
SNS junctions. Indeed the resistance and geometrical 
capacitance of SNS junctions (R a few Ohms [22|, and C 
roughly 10~ 16 fF) would correspond to Q -C 1. In this 
so called overdamped regime, the expected I(V) depen- 
dence is V = \Jl 2 — 1% , and therefore no sharp switching 
or hysteresis should be observed [TTI | . Experimentally we 
find I s / I r ~ Q between 3 and 7, see table I, and in other 



groups hysteresis is always observed at low temperature 
(Dubos et al. find 8, Anthore finds between 5 and 6 [HI). 
A similar discrepancy has already been pointed out in 
superconducting weak links by Song [23| . The large hys- 
teresis was explained by the fact that the pair relaxation 
time in the weak link was no longer given by RC but by 
h/A. Following this analysis, we replace the RC time 
by the diffusion time of the Andreev pairs in the normal 
region tr> — h/E t h- Comparison with the RSJ model is 
then achieved by replacing the capacitance with an effec- 
tive capacitance C e // = td/R. This yields an expression 
for the quality factor: Q = I c /I r = \j2e/hRfqI c /TB = 

\j2eRNl c / Etii = \/2b. Table U shows that the quality 
factor is indeed of the order of this value. 





SIS 


SNS 




A 


bE Th 


Ej 


(ft/e)Je 


(h/e)I c = ±^E Th 


'RC' 


RC 


h/E Th 


E c 


e 2 /C 


27v£-E Th 


hujp 


^/h2eI c /C 


V2bE Th 


Q 


ujpRC 


V2b 


k B T* 


hui p /(2n) 


V2b pi 



TABLE II: Comparison between quantities determining the 
switching process, for SIS and SNS junctions. Here Rq = 
h/e 2 and b = cRnIc/ 'E Th . 

Note that heating of the sample once the SNS junc- 
tion has switched to the dissipative state will also cause 
hysteretic I V curves: The power injected in the nor- 
mal wire may not be entirely dissipated by the substrate 
and superconducting contact phonons at low tempera- 
ture, and heat conduction by the superconducting con- 
tacts is weak because of the superconducting gap. It is 
therefore probable that the electron temperature as the 
current is ramped down after switching is higher than be- 
fore switching, so that the experimental retrapping cur- 
rent may be smaller than the intrinsic retrapping current. 

Since both heating and intrinsic charge relaxation de- 
pend on the wire length (albeit with different power 
laws), it is difficult at this stage to determine exactly 
the relative contribution of each effect to the hysteresis. 
However the very shape of the IV curve is a clear indica- 
tion that these SNS junctions are intrinsically hysteretic, 
since a non-hysteretic (overdamped) junction would have 
a square-root dependence V = R\jl 2 — I 2 of the V(I) 
curve [ll| . in contrast with the sharp jump observed at 
the transition in our junctions. Furthermore, measure- 
ments of the noise spectrum of long diffusive SNS junc- 
tions by Hoffmann et al. [2(| suggest that the electron 
temperature is less than 500 mK, even in the retrapping 
branch. 

It is therefore safe to conclude that although electron 
heating may contribute to the hysteresis in the IV curve, 
intrinsic hysteresis is also present and most likely due to 
the relatively slow dynamics of the phase and charge in 
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the long normal wire. 

The inset of fig. [5] shows that the hysteresis disappears 
at higher temperature. This could be due to a lower I c 
and to a more efficient heat dissipation at higher tem- 
perature, and therefore decreased heating. But it may 
also be explained by the fact that the quality factor Q 
depends on temperature via the square root dependence 
of the plasma frequency oj p upon I c (see table II) . In the 
case of this sample Nb750, the switching current is di- 
vided by 8 between low temperature and 1 K, and the 
quality factor should thus loose a factor y/8 ~ 3, bring- 
ing the quality factor to 1 at IK, corresponding to strong 
damping and a non hysteretic curves with square root 
V(I) dependence, as is observed (see Fig [5]). 



b. Fluctuations of switching current The dynamics 
of the SNS junction also determine the probability dis- 
tribution of the current values at which the junction 
switches from its superconducting to normal state. We 
have measured histograms of the switching current, fol- 
lowing the technique used by Devoret et al. 30]. The 
current is repeatedly swept (at a frequency of / = 65 Hz) 
from a small negative current to / ma x > 7 C , and the cur- 
rent at which a finite voltage appears is recorded. Fig[7] 
shows a histogram of the switching current for a Nb/Au 
ring at §/$ = (2n + 1)tt. The histograms of Al/Au 
junctions are similar. As expected from a switching phe- 
nomenon, the histogram is asymmetric, falling of more 
sharply at the higher currents. The histogram is narrow, 
with SI S /I S roughly 1 per thousand. 



Following the previous description of the hysteresis, 
in which the geometrical capacitance C is replaced by 
an effective one td / R, we can extrapolate the SIS the- 
ory to the SNS case. All the characteristic scales for 
SIS junctions and their correspondence for SNS junc- 
tions are listed in table |TTJ Garg [2^| has calculated 
the width of the switching histogram of SIS junctions. 
This width depends on the damping and the origin of 
the escape rate (thermal or quantum). In the under- 
damped regime, which is our situation (see table), at 
high temperature the particle should escape by ther- 
mal activation and the width is given by SI S /I S ~ 
(T/Ej) 2 / 3 but below a temperature T* = fiujp/2-nkB, 
of the order of Etk-, the particle escapes by quantum 
tunneling and the width of the histogram is given by 
5I S /I S pa (E T h/ Ej) 4 "/ 5 . Extrapolating these results to 
the SNS junction, the expected width is proportional 
to {T/Ejf/ 3 = {2TTR N /R Q .T/E T h) 2/3 in the thermally 
activated regime, and proportional to (Erh/Ej) 4 / 5 = 
{2-kRn I 'Rq) a / 5 in the quantum regime. The predicted 
widths are therefore of the order of 6I S /I S = 1 per thou- 
sand, as measured experimentally at low temperature. 
The 2/3 power law however in the thermally activated 
regime is not really found experimentally (see Fig. [?])■ 
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FIG. 7: Left pannel: Switching current 7 S and width of the 
histogram (Standard deviation SD) as a function of tem- 
perature for the short Nb dc squid (Nbsq750), at $ = 0. 
Right pannel: ratio SD/I S , compared to a linear and T 2//3 
dependence. Inset: Histogram of the switching current for 
$/$o = (2n + 1)tt for the shortest junction Nb/Au ring. The 
total number of counts is 2000. 



E. Temperature dependence, Thouless energy and 
normal state resistance 



The critical current is predicted [12|, [l3| to reach 
eR N I c {T = 0) = 10.82 E Th in the limit of infinitely 
long junctions (A/Exh — oo), and to be given by 
eR N I c (T = 0) = bE Th , with b smaller than 10.82 but 
within a few percent as soon as A/Eth > 10. The ra- 
tios b — cRj^Isj Etk at the lowest temperature (16 mK 
for Nbsq and Alsq samples, 60 mK for Alw samples) are 
listed in table [T] and are within 50% of this value. The 
main uncertainty in calculating the ratio however is in 
determining the normal state resistance. The Thouless 
energy Etk is extracted with greater confidence from the 
temperature dependence of the switching current, as ex- 
plained in the next paragraph. 

The switching current can thus be defined (rather 
equivalently) as either the average or the most proba- 
ble value out of a series of measurements, and is plotted 
for 2000 measurements as a function of temperature in 
Fig. [5] As is clear in the figure, the decay of switch- 
ing current above a few hundred milliKelvin is roughly 
exponential: I s ss e~ T / aETh , w jth a of order 9 [l3|]. We 
used an approximate high temperature expansion (Equa- 
tion (1) of Dubos et al. [l6j |. valid when fc^T > 5Exh) 
to determine the Thouless energy and the normal state 
resistance (dashed lines). The full temperature depen- 
dence, obtained by numerical integration of the Usadel 
equations using the determined Thouless Energy, is also 
plotted as a continuous line of the figure. For both Al and 
Nb samples the high temperature decay of the switching 
current is well described by the Usadel calculation. At 
low temperature however, the predicted switching cur- 
rent is larger than found experimentally (by roughly a 
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factor two, for instance 100 fiA instead of 70 \iA for the 
longests Nb dc squid (Nbsql200), and 500 \i A instead of 
250 [i A for the shortest dcsquid (Nbsq750)). The devia- 
tions from the theoretical prediction occur below 400 mK 
for the longest junction, and below 200 mK for the short- 
est one. The switching current saturates at lower tem- 
peratures. We have no explanation for the deviations 
from theory at temperatures which are relatively high: 
although the electronic temperature may be higher than 
the phonon temperature of the mixing chamber at the 
lowest temperatures (16 to 50 or 100 mK), it is improba- 
ble that the electron temperature be very different than 
the phonon temperature at temperatures higher than 100 
mK [3. 

The fitting parameters Rn and Eth are given in Table 
1, as are the value of the ratio eR^Iaj Etk- The nor- 
mal state resistance Rn extracted from the fit is some- 
what (50%) smaller than that extracted from the dV/dl 
curves, and the eRwIs/ Etk product is also roughly 50% 
less than maximum value 10.82, expected at zero tem- 
perature, for an infinitely long normal wire and a perfect 
transparency. 

The role of non perfect transparency of the SN inter- 
face has been investigated in ref [l5j. Numerical integra- 
tion of the Usadel equations shows that the main effect of 
a non transparent interface is to renormalize the Thouless 
energy and to decrease the region of current saturation: 
thus an imperfect interface could not explain the slower 
variations with temperature of the switching current at 
low temperature. 



F. Magnetic field dependence 

The variations of the switching current in a perpen- 
dicular magnetic field is central to our investigation. As 
for conventional SIS dc squids, the SNS dc squid samples 
have a periodically oscillating switching current, whose 
period corresponds to a flux quantum <j> = h/2e through 
the loop area, as shown in figures I9I10I11I12I and [13] 
There are however features of the field dependence that 
have not been reported before, for instance the absence 
of the Fraunhoffer pattern in the decrease of switching 
current at larger field: the decay is monotonous, prac- 
tically gaussian with a field scale of one flux quantum 
through the normal wire. 
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FIG. 8: Switching and retrapping current as a function of 
temperature for both Nb dc squids, at maximal switching 
current, comparison between experiment and theory. Filled 
dots and squares are the data points. Red dashed lines are 
fits to the high temperature formula for the proximity effect 
(see Eq (1) of Dubos et al., from which the Thouless energy is 
deduced. Black lines are the numerical solutions to the Usadel 
equations using the Thouless energy deduced from the high 
temperature fit. The hysteresis disappears (i.e. the switching 
and retrapping current are equal) at a temperature close to 
llE Th . 




FIG. 9: Switching current of Alsql900 SNS dc squid as a 
function of magnetic field at T= 16 mK (circles), and fit to 
Equation [3] (line) showing the dc squid-like behavior. Here 
h = 2.1 iiA and h = 2.7 pA. 



1. Modulation of the switching current by the magnetic 
field: the DC SNSquid 



We first focus on the proximity SQUID interference 
patterns, presented in figures [9] and 1 101 and compare 
them to the theory of a regular DC SQUID made of two 
tunnel SIS junctions in parallel 11!]. The critical cur- 
rent I c , which is the switching current in the absence of 



fluctuations, is given by: 



It = (h - hf + 4/i/2 cos 2 (^^ 

<£>() 

= If + /| + 2/1/2 cos(2tt-5-) 



(2) 
(3) 



8 




FIG. 10: Dots: Switching current over three magnetic field 
periods, for the longest Nb dc squid Nbsql200, measured at 60 
mK. Lines are the simulated curves corresponding to the same 
maximum switching current of 70 fiA: black dotted curve, a 
non inductive squid with identical junctions (same critical 
current); red dash-dotted curve, a non inductive squid with 
different junctions; and green line, a squid with identical junc- 
tions but two inductive branches, with identical inductance 
L=6 pH. 

with ii and I2 the critical currents of the 2 junctions, and 
$ the magnetic flux threading the SQUID loop. This crit- 
ical current thus oscillates periodically with field with a 
period of one flux quantum <!>o = h/2e through the squid 
loop. The oscillation amplitude is maximum in symmet- 
ric squids, and is reduced in squids whose junctions have 
different critical currents. The relative amplitude in non 
inductive squids is given by 

jrnax _ jmin ^ 2min(/ 1 ,/ 2 ) 

Vas ~ I™ ax ~ h + h 

Such a dc squid behavior is indeed demonstrated in our 
SNS dc squids, as shown in Figs. H [TOJ HQ H ESI Super- 
currents up to 330 \iA are modulated by 45 to 88%, see 
table 1. The incomplete modulation, which corresponds 
to a ratio I1/I2 between 1.3 and 3.4, can be attributed 
to slight differences in the geometry of the two normal 
wires, to which the critical current is very sensitive. In- 
deed, the critical current of a long SNS junction scales as 
the inverse cube of the normal wire length. 

The inductances L\ and L 2 of the 2 branches of the 
SQUID also come into play in determining the shape of 
the modulation curve, since they act to screen the exter- 
nal flux $ ex t by <5$ = L±ii — £2*2 where i\ and 12 are 
the currents through the 2 branches of the SQUID (and 
ii + *2 = I)- In symmetric squids with L\ = L 2 , the 
screening induced is zero for fields corresponding to the 
maximum switching current, since the flux induced by 
the screening current in one branch is exactly canceled by 
the flux induced by the same current in the other branch. 
The maximum effect of the screening flux produced by 
the inductance corresponds to half a flux quantum in the 
loop &ext/&o = 1/2 mod 1. This explains deviations 
from the simple cosine dependence seen in some samples 
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FIG. 11: Sample Alsq900, comparison between the positive 
switching current in the up sweep (increasing current, red 
square symbols) and the negative switching current in the 
down sweep (decreasing current, black round sumbols). The 
black curve with small open dots is the transformation of 
the down switching curve, via the time reversal operation 
T(lf n (H)) = —lf n (—H), and is seen to be superimposed on 
Ig P (H), as expected from the global time reversal symmetry 
of the experiment. 

(see for instance Fig. [TO]). In some cases, corresponding 
to L\ ^ L2, asymmetries with respect to reversal of mag- 
netic field or current are clearly visible (see Fig. [TTj) . The 
switching currents I" p and lf n measured on the positive 
and negative current branches of V(I) (while increasing 
the absolute value of current in both cases) differ at finite 
magnetic field: 

I?(-H) jt I7{H) (5) 
If n {-H) + lf»(H). 

However, as expected from the global time reversal 
symmetry of the experiment, I^ P (H) = —lf n (—H). 
These asymmetries have previously been calculated in 
[3l|, and measured in S/I/S dcSQUIDS [Tp . They were 
also exploited by Clarke and Paterson Q in larger area 
SNS dc squids: the authors controllably adjusted the 
asymmetry of the two superconducting branches and ver- 
ified the corresponding asymmetric I s (&) curves. The 
current modulation in those SNS dc squids was less than 
5 percent. As expected for self inductance effects, the 
asymmetry and anharmonic distortion of the periodic os- 
cillations are largest in samples with the largest critical 
current. 

One also expects the intrinsic non harmonicity of the 
current-phase relation of an SNS junction to cause a dis- 
tortion in the switching current Vs field curve of an SNS 
dc squid: the amplitude of the m-th harmonic in the (f>o 
periodic current phase relation is predicted to vary like 
1/m 2 pjj and should therefore be detected at least up 
to the 3rd order. However as is shown in the Appendix, 
the critical current of a dc SNSQUID with 2 symmetri- 
cal junctions is insensitive to even harmonics, so that the 
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FIG. 12: Switching current as a function of magnetic field for 
Al SNS dc squid Alsql500, at 17 mK. The line is a Gaussian. 

harmonics content of the current-phase relation does not 
easily show up. To specifically test the anharmonicity 
of the current-phase relation in a long SNS junction, a 
strongly asymmetric setup with one reference SIS har- 
monic junction like the quantronium 18], would in prin- 
ciple be more adequate. However, it would still be sen- 
sitive to self inductance effects for large supercurrents. 
Therefore the best way to measure the various harmon- 
ics of the current phase relation of long metallic junctions 
is to directly measure the magnetic orbital response of a 
simple SNS ring, as was done recently using a magnetic 
Hall probe technique 



2. Decrease of the critical current at high field 

We now discuss the decrease of the switching current 
at large field scales. In single junction devices, the extinc- 
tion of the switching current with field is straightforward. 
In the dcSNSquid samples, whose switching current os- 
cillates on a small field scale, we consider the extinction 
at large field of the envelope of the oscillations. Where 
necessary, we follow the behavior with field of the max- 
ima of the oscillations. We will in the following neglect 
the effect of the asymmetries in the squids, and consider 
that the behavior in field of these maxima is equivalent 
to the behavior in field of a single SNS junction. 

The first observation is that we do not observe the 
Fraunhoffer diffraction pattern encountered in many ex- 
periments on wider and shorter SIS and SNS junctions 
0; E3, EH- Instead, we find a monotonous decay with 
a typical field scale of a flux quantum through the nor- 
mal metal wire. For most of our samples, the decay is 
approximately gaussian. 

Qualitatively, in a semi-classical picture, the decay can 
be seen as due to dephasing induced by the magnetic 
field: different diffusive Andreev paths in the normal 




FIG. 13: Switching current of the short junction Nb/Au dc 
squid Nbsq750 as a function of magnetic field, at 60 mK. 
Black line is a Gaussian. 

metal acquire different phases. If the junction is wide, 
the phase along the NS boundary also varies, but in the 
limit of a normal wire much longer than it is wide, only 
the phase acquired during the propagation across the nor- 
mal metal is relevant. 

Quantitatively, the average of the phase factor 0(C) = 
J c Adl I $o acquired along a path C has been determined 
analytically in the ID limit (normal wire length L much 
greater than width W), using the solution of the field 
dependent diffusion equation [42]. This leads to 

TT H 

i c {H)=m ff° H (6) 

Slim V3 H 

where Hq = <&q/WL corresponds to a flux quantum 
through the normal part of the SNS junction (and <J>o = 
h/(2e)). 

This calculation does not take into account the modifi- 
cation of the density of states in the normal wire. These 
experiments have also stimulated a full treatment of the 
proximity effect via the numerical resolution of the Us- 
adel equations, for any ratio of the wire length versus 
width. The magnetic field enters in these equations via 
a depairing time th — 6h 2 e 2 / (e 2 Dw 2 H 2 ) . Cuevas et al. 
[TH have found that in the ID limit (L > W), the so- 
lution can be fit by a gaussian dependence at low fiels 
(H/H < 3). 

I C (H) = J c (0)exp- [a(H/H ) 2 ] . (7) 

the coefficient a is worth 0.24 at zero temperature, and 
decreases as T increases. 

For low fields (H/Hq < 2), the semiclassical calcula- 
tion equation [5] is also well approximated by a Gaussian, 
but the coefficient in the exponent is 0.55 instead of 0.24. 
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FIG. 14: Field dependence of the normalized switching cur- 
rent for all samples at lowest temperatures. For the samples 
with the squid geometry, the filled data points consist of the 
local maxima of the full (modulated) curves. Also plotted 
with open symbols are the theoretical predictions of the Us- 
adel equations at zero temperature, for four of the samples 
taking into account the exact aspect ratio of the junction. 
The prediction of the ID semiclassical calculation Eq. [6] is 
the continuous line. The vertical scale is logarithmic, and 
the horizontal axis is the square of magnetic field normalized 
by the field Ho corresponding to a flux quantum through the 
wire surface. 



In SNS junctions with wider wires (W > L), Cuevas 
et al. [15] find that the critical current does not decay 
monotonously with field but that oscillations appear, and 
that for very large junctions the well known Fraunhoffer 
pattern is recovered. 

The ratio L/W of our samples is respectively 2 and 3 
for the short and long Nb samples, and exceeds 7 for all Al 
samples (see table|TJ . To compare experiment and theory, 
we thus computed the predicted field dependence of the 
maximal switching current, using the Usadel equations 
with the experimental aspect ratio of the junctions. The 
curves for both Nb SNS dc squids and two Al SNS sam- 
ples are plotted in Fig. [TU It can be seen that the decay 
of the switching current is only qualitatively described 
by theory. In particular theory predicts that the widest 
Nb/Au/Nb junction should exhibit relatively large os- 
cillations of the maximum switching current with field, 
which is not found in the experiment. 

One possible explanation for these discrepancies is the 
screening of the magnetic field by the superconducting 
contacts, as well as by the normal parts. Screening by the 
superconductor focuses the magnetic flux lines through 
the normal, leading to a smaller Ho. Screening by the 
normal metal on the other hand tends toward flux ex- 
pulsion and a larger Hq. We have not yet included these 
effects in the calculations. 

Fig. [Pol shows the field dependence of sample Alsql500 
at several temperatures, the temperature dependence of 
the Gaussian coefficient, and the comparison with the 



FIG. 15: Field dependence of the switching current, as a func- 
tion of temperature. Left panekMaximum switching current 
(extracted from the modulated curves) as a function of mag- 
netic field for Al SNS dc squid Alsql500, at temperatures of 
(from top to bottom) 16, 55, 110, 168, 233 and 304 mK (sym- 
bols). The vertical scale is logarithmic, and the horizontal 
axis is the square of the reduced magnetic field. Dashed lines 
correspond to a Gaussian decay. Right panel: Full symbols, 
temperature dependence of the coefficient of the Gaussian de- 
cay (coefficient a of Eq. extracted from the data, compared 
to the prediction of the Usadel equations (empty symbols). 



coefficients predicted by the Usadel equations. Here also 
experiment and theory differ, a fact that may be due to 
the temperature dependence of screening of the magnetic 
field. 



3. Reentrant behavior at low field and low temperature 

Certain samples exhibit a surprising increase of the 
switching current at low magnetic field (below 20G, see 
fig fTGf and at temperatures below 150 mK. A qualita- 
tively similar effect but on a very different field scale, was 
already observed on narrow superconducting wires made 
of MoGe and Ge [34[ and explained by phase breaking of 
Cooper pairs induced by spin flip scattering on magnetic 
impurities. Such spin flip scattering disappears when the 
magnetic moments are polarized under a magnetic field 
such that hbB = ksT. This would correspond to a field 
of a few hundred Gauss in the present experiment, more 
than an order of magnitude larger than observed. Similar 
anomalies at low field were also reported by Dynes et al. 
[35l ] in narrow supercondcuting wires, and attributed to 
phase fluctuations of the superconducting order parame- 
ter. Other possible explanations include: 

- Flux trapping in the superconducting wires during 
cooling of the samples. This can probably be excluded 
since no irreversibility or hysteresis is associated to this 
unusual field dependence. 

- Quantum interference in the normal part of the junc- 
tion. They are known to lead to weak localisation cor- 
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FIG. 16: Switching current as a function of magnetic field 
for the three Al/Au SNS wires, at 17 mK. Switching current 
of wire Alw900 has been divided by 2 for better comparison 
with the two other samples. Two of the three wires show a 
reentrant behavior. 

rections of the conductance in the absence of proximity 
effect. The correction is of the order of e 2 /h, which for 
a wire of roughly 1 f2 _1 corresponds to ^ ~ 10~ 4 , i.e. 
much less than the 10 or 20% observed in the experiment. 
Moreover, the reentrance disappears around the tempera- 
ture corresponding to the Thouless energy, to which weak 
localisation is not sensitive. 

- A last possibility could be a magnetic orbital effect, in 
the form of quasiparticle persistent current loops, giving 
rise to strong paramagnetic orbital contributions similar 
to what was observed in hybrid N/S cylinders at very low 
temperatures [36l l37l|. 

It is noteworthy that this reentrance of I S (H) is only 
present in the aluminum S/N/S junctions or SQUIDs, 
which have an abrupt interface between the normal and 



superconducting regions, and was never observed in the 
Nb junctions, in which the Au is also present under the 
Nb electrodes. It also occurs only in normal junctions 
whose length is between 1.2 and 1.5 fim, and was not 
detected in shorter or longer wires. 



G. Conclusion 

The SNS dc squid resembles its SIS counterpart, with 
a strong modulation of its switching current by the mag- 
netic field, with a periodicity of one flux quantum in the 
loop. But it also has unique features, such as a large 
switching current due to the low impedance of the normal 
metal. The IV curves and switching current histograms, 
because they depend on the intrinsic response functions 
of the normal metal wire, suggest the importance of in- 
vestigating the dynamics of the Andreev states. In con- 
trast to the SIS dc squid, the length and width of the 
normal part of the junctions can be varied at will. In the 
samples presented in this paper, the aspect ratio of the 
normal wires places them in the ID or narrow 2D limit. 
We have shown that in this limit the decay of switching 
current at high field does not show an oscillating inter- 
ference pattern typical of wider junctions, but instead 
decays monotonously in a Gaussianlike manner. Finally, 
some samples have a surprising increase of switching cur- 
rent at low field, which is still not well understood. 
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I. Appendix: Appearance of higher harmonics in a 
dc SNSquid 

The current phase relation in a SNS junction is not 
simply a sinusoidal relation but is predicted [l2j to con- 
tain higher harmonics: 



J S ivsW = X] I " sin (^) ; (8) 

n=l 

with in the case of a long junction 

{-\r^eE Th 
c ~ Ejv(2n + l)(2n-l)' 1 ' 

The supercurrent through an SNS dc squid with two 
identical SNS junctions in parallel thus verifies: 

j t = (E 7 " s[nn ^ + 0) + E 7 " sin n ^ - *))) ( 10 ) 

n n 

= 2 E I™ sm n ^o cos n<p) , (11) 

n 

where 4> = e&/h and So + 4> and So — <fi are the phase 
differences at each of the junction. Around So = (2p+l)-| 
(which is no longer the position of the maximum, but is 
close to it), it is easy to see that only the odd harmonics 
are non zero, so that the critical current is given by 

jSns = 2 J- 7«|cos^|. (12) 



